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Abstract 



We consider scalar field perturbations about asymptotically Lifshitz black holes with dynamical 

^^' exponent z in D dimensions. We show that, for suitable boundary conditions, these Lifshitz black 

r^ ' holes are stable under scalar field perturbations. For z = 2, we explicitly compute the quasinormal 

mode frequencies, which result to be purely imaginary, and then obtain the damping-off of the 

Cn . scalar field perturbation in these backgrounds. The general analysis includes, in particular, the 

^^ \ z = 3 black hole solution of three-dimensional massive gravity. 
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'^ ■ 1 Introduction 

m 

^ ! 1.1 Lifshitz holography 

It has been conjectured in [1] that the non-relativistic Lifshitz-like fixed points admit a dual descrip- 

yj . tion in terms of holographic gravity models, mimicking what happens in the standard AdS/CFT 

Vh ■ correspondence for the case of relativistic conformal field theories [2]. The proposal in [1] is part of a 

■ ■ ■ more general program, whose goal is to generalize AdS/CFT correspondence to the area of condensed 

matter physics. The importance of having a tool to describe strongly-correlated condensed matter 

systems is undoubtable, and this is the reason why much attention has been focused on this area in 

the last four years. 

In this context, a Lifshitz fixed point refers to a non-relativistic system described by a field theory 
exhibiting anisotropic scale invariance under the transformation 

x^\^ X , t^\^^ t (1) 

with z G M>i in the spacetime coordinates, but no Galilean invariance. Systems with such a be- 
havior appear, for instance, in the description of strongly correlated electrons. These systems are 
characterized by the value of the dynamical exponent z that parameterizes the scale transformation 
p!]). The proposal in [1] is that, analogously to the case of relativistic conformal field theory, the 



{D — l)-diniensional Lifshitz fixed points would also admit a geometric dual description in D dimen- 
sions. The geometric realization is such that the D-dimensional asymptotically Anti-de Sitter spaces 
of the standard AdSD/CFT£)_i picture get replaced by spaces whose metrics asymptote the geometry 

ds^ = --p^dt^ + -^dr^ + r^ dxl,_o (2) 

at large r, where dx'jj_2 is the metric of the {D — 2)-dimensional plane M , I is a positive parameter 
with dimensions of length, and z is a real parameter that ultimately coincides with the dynamical 
exponent of the scale invariant theory that one is trying to holographically describe. Then, spaces ([2]), 
dubbed "Lifshitz spaces", happen to realize geometrically the scaling symmetry ([1]), which is certainly 
a symmetry of the metric above provided one also rescales the radial coordinate as r — )• X~^ r. In this 
generalized framework, the standard AdS£)/CFT£)_i correspondence is recovered in the particular 
case z = 1, for which ([2|) coincides with AdS^) in Poincare coordinates. This suggests a natural way 
of prescribing the holographic dictionary by mimicking the recipe of [3l |3j and adapting it to generic 
values of z. This holographic correspondence was thoughtfully studied, for instance, in references 
[H [71 [51 [HI [To]. As seen in [1] and in its sequels, the fact of having -z / 1 introduces new features 
in the boundary field theory correlators. For instance, the fact of having replaced AdS space by the 
Lifshitz space introduces a dependence of z in the formula that gives the scaling dimension of the 
operators in the dual theory: In the case of a scalar field of mass m in a D-dimensional Lifshitz space 
of dynamical exponent z, the scaling dimension A of the dual operator is given by the formula 

A"^ -{D + z-2)A- rr?l^ = 0, (3) 

and its value corresponds to the asymptotic behavior of the bulk field at large distance, namely 

$(x, t, r) ~ r-^ ^q{x, t) + O(r-^-i), (4) 

with the boundary configuration ^Q^x^t) acting as a source for a dual field from the boundary theory 
point of view. Finiteness of the scalar field action on ([2|) demands 

2A> D + z-2, (5) 

that is, it demands the following z-dependent generalization of the Breitenlohner-Freedman bound 

vr? > -'- ^^^^. (6) 

1.2 Correlation functions 

Correlation functions in the dual theory can then be computed by following the standard holographic 
prescription, which requires to calculate the on-shell regularized Euclidean action with appropriate 
boundary asymptotic. To do this, one first performs a Wick rotation, defining r = it. It is also 



convenient to define inverse radial coordinate as u = 1/r, for wliich tlie boundary turns out to be 
located at u = 0. Tfien, one solves the scalar wave equation on the Euclidean background, whose 
solution, in the near boundary region, namely for small values of n, behaves like 

$(x,r,n)~u^+ ${,+^(3;,t)+u^- $[,"^(x,r), (7) 

where A-t are the solutions to (l3|), cf. (JH). 

To calculate the boundary two-point function one needs the bulk-boundary propagator. One writes 
the bulk field configuration in terms of the Green function G{x, r, u] x', r', 0); namely 

*(...,u)=/..'..'.(/,/.0)G(...,u;.',.'.0). (8) 

In the momentum space the convolution simply reads 

^{k,uj,u) = G{k,oj,u)'^{k,uj,id). (9) 

Then, one evaluates the bulk action, first integrating by parts and introducing an infrared cut-off. 
Namely, 

//•oo /• 

In the momentum space, one finds 

5[$] = I dkduj{^{k,uj,f))F{k,u) $(-A:,-w,0)) (11) 

with the flux function F{k^ u) being defined as follows 

r _ ~ -1 oo 

F{k,u:)= G{-k,-io,u)^g^^duG{k,u;,u) . (12) 

J e 

This eventually gives the boundary two-point function in the momentum space 

F{k,io) = {0{k,uj)0{-k,-uj)) (13) 

where only the flux factor evaluated on e — t- contributes, as the propagator G{k,uj,u) vanishes when 
u — )• oo. Nevertheless, the quantity one is interested in is the two-point function in position space, 
which is obtained by Fourier transforming the expression above. This requires a careful treatment of 
the contact terms when Fourier transforming the distributions [11] and taking the large separation 
limit. 

In the particular case of strongly coupled z = 2 Lifshitz points in D — 1 = 3 dimensions at zero 



temperature, it was shown in [T] that the two-point function in the position space exhibits the form 

(0(X1, tl)0(x2, t2)) ^ I ^—T^A + - > (14) 

where c is a non-vanishing constant, 0{x,t) is an operator of the boundary theory, dual to the bulk 
operator $(x,t,r), and where the ellipses stand for contributions that are subleading in the large 
separation limit and contact terms. Notice that, in the case D = 4 with z = 2 and tti = 0, we have 
A = 4. 

A remarkable feature of (|14p is that it does not involve further suppressions coming from overall 
functions of the scale invariant quantity |xi — X2p/|ti — ^2! as one could have expected. As discussed 
in [1], such a suppression could have resulted in ultralocal behavior of the two-point function. Taking 
a closer look at the computation of (|14p one observes that, in general, the reason why no ultralocal 
behavior is found in the boundary correlation function, is the existence of non-analytic contributions at 
small momenta. In the position space, such non-analytic contributions yield a non- vanishing expression 
when \xi — X2I > 0. If only positive integer powers of the momenta were present, then the Fourier 
transformation would give only contact terms. 

1.3 Finite temperature and the quasinormal modes 

Finite temperature regime, on the other hand, is also modified with respect to the standard z = 1 
case. We know from AdS/CFT that adding finite temperature in the boundary conformal field theory 
corresponds to considering a black hole in AdS space. Similarly, in the case of anisotropic scale 
invariant theories, adding finite temperature would correspond to considering in the bulk a black hole 
solution that asymptotes to ([2]) at large r. 

The two-point function at finite temperature is computed along the same lines as in the zero- 
temperature case: Basically, one solves the scalar field equation in the Euclidean black hole background 
imposing incoming-modes boundary conditions at the horizon. Typically, the finite-temperature two- 
point function would exhibit an exponential damping-off produced by temperature effects. For in- 
stance, in Z) = 4 and z = 2 this was explicitly computed in [12], and the result was shown to have the 
form 

— \/47rX|x'i — xol 

(0(X1, ii)0(x2, t2)> ^ -. ^7^ + ... (15) 

Fi - a;2r/^ 

Apart from being necessary to solve the field equation on the black hole background to compute 
the two-point function, the black hole quasinormal mode frequencies also have a direct interpretation 
from the boundary point of view: Considering a large black hole in the bulk corresponds to considering 
a thermal state in the dual boundary theory, and the decay of the scalar field in the bulk is directly 
associated to the decay of a perturbation of such a thermal state, whose relaxation time is given by 
the imaginary part of the quasinormal frequency u}{k), |14[ l22j. 

For the case of z = 2 Lifshitz black holes, we will explicitly solve the quasinormal modes, analyti- 
cally computing the frequencies and then obtaining the damping-off of the scalar field perturbations. 
Then, adapting the computation that holds for asymptotically AdS black holes to the case of generic 



z, we will also show that the imaginary part of the frequencies for the scalar field perturbation about 
the Lifshitz black hole backgrounds is always negative, proving the stability of these geometries under 
this kind of probe perturbations. 

2 Lifshitz black holes 

In this section we will construct a higher-curvature gravity theory that admits z = 2 Lifshitz black 
holes in arbitrary dimension D. We begin by considering the following ansatz 

^'D = -p^F (0 dt^ + ^:2^^ + r'dxl_^ , (16) 

asking the metric function to obey F [r ^ +oo) — t- 1 and have a simple root for r > 0. 

Analytic Lifshitz black hole solutions are not easy to be found, and, consequently, constructing 
finite temperature gravity duals has so far required the introduction of strange matter content whose 
physical interpretation and theoretical motivation are unclear. An alternative way of finding such 
a Lifshitz black hole solution is considering carefully-tuned higher-curvature modifications to the 
Einstein-Hilbert gravity action. This has been done, for instance, in references [I6]-[2T]. Here, since 
we will be concerned with the Lifshitz black hole solutions with z = 2, we will present the simplest 
higher-curvature gravity Lagrangians that admit asymptotically Lifshitz black holes with z = 2 in 
arbitrary dimension D. As we will see, this amounts to considering terms in the gravity Lagrangian 
that are of cubic order in the Riemann tensor. 

Then, we look for higher-curvature Lagrangians that admit metric (|16p as solutions with z ^ 1. 
Since the most conservative generalization of Einstein gravity in higher dimensions. Lovelock theories, 
excludes static Lifshitz black holes [23], it is necessary to explore theories of gravity with higher order 
field equationq^. In such a setup, black hole solutions have recently been reported; in particular, in the 
quadratic D = ?> massive gravity of |13j a black hole with z = ?> and F {r) = 1 — r^jr^ was found in 
reference |16j . In reference |18] the existence of such solutions in generic quadratic gravity theories in 
D dimensions was explored. However, in the setup of [18] it is not possible to achieve the value z = 2 
for the dynamical exponent, and so a natural next step is to include cubic terms in the Lagrangian. 

We are interested in a theory that would result valid in D = 3 as well, so that we consider only 
higher-curvature terms constructed with the Ricci tensor R^y. In three dimensions, the Weyl tensor 
is zero and the Riemann tensor is fully determined in terms of its trace. Then, we only have two 
quadratic invariants, namely 7^?-^n = i?^ and T^^g) = R^uR'^n-, and three cubic invariants, 7^Q^ = i?^, 
7^?2) = RR^uR^"^, and T^%) = R'^uR'^pR^'fi- For D > 7, there are eight linearly independent, cubic 
scalars which can be constructed out of three Riemann tensors without involving its derivatives (see 



^The black holes and wormholes in vacuum, found in [T^ and [20] respectively, exist only for fixed values of the 
couplings of Lovelock gravity, at which Birkhoff 's theorem are circumvented. 



e.g. |24j ) . Here, for simplicity, we consider the simplest model defined by the action 

1=1 d^xV=^ iaR-2A + Yl ^i'K) + iZ ^*^0) I (1^) 

Already in D = 3 one observes that this Lagrangian admits Lifshitz black hole solutions with 
z = 2. More generally, one verifies that this Lagrangian happens to admit asymptotically Lifshitz 
black holes of the form (J16p with z = 2 for all values of D. Notice also that, since the action contains 
only traces of the Riemann curvature, this theory admits all Einstein manifolds as solutions, for which 
R^u = ^g^iu with A being fixed in terms of the couplings. It is not hard to show that the field equations 
coming from (|17p admit the following black hole metric as a solution 

dsl = --[l-^]dt' + — ^ + r^dxl_^ , (18) 



l^\ r^J (r2 - rl) 

provided the coupling constants A, 6j, and q obey specific relations that we write in the Appendix. 
Here, the inverse of the Newton constant a and the parameter / can be regarded as two arbitrary 
parameters. Theory (J17p also admits Lifshitz black hole solutions with other values of z. 

Black holes (jlSp present a Hawking temperature T = r'^/{2irl'^). In the metric, r_|_ is an integration 
constant which represents the location of the event horizon, r = r_|_. This integration constant is an 
actual physical parameter provided at least one of the flat directions is compactified, which results in 
breaking the scaling symmetry ([1]). On the contrary, if all the flat directions are non-compact, then 
r_|_ can be settled to 1 by a simple change of coordinates. 

Having given a particular example of a theory that admits z = 2 Lifshitz black holes in arbitrary 
dimensions, we can proceed and study the scalar response on these backgrounds. 

3 Quasinormal modes 

3.1 Quasinormal modes for z = 2 black holes 

As pointed out in [12], the scalar field response in backgrounds (jlSp turn out to be exactly solvable 
in terms of hyper geometric functions. This is relevant, for instance, to compute boundary correlation 
functions. This is precisely the computation done in reference |12j . where the finite temperature 
two-point function was shown not to exhibit ultralocal behavior either. Here, we are interested in a 
different problem, namely that of computing quasinormal mode frequencies on background (jlSp . To 
solve this, we consider a massive scalar field on the black hole metric psp 

,2\ ^ _ ^ o / / — -„uua /T,\ _2; 



{n-m^)^ = ^=^f,{^/^g^"'^^^)-m^'^ = , (19) 

and assume separability, namely 

$(f,t,r)=e-*"'*+*^-^i?(r) , (20) 



where k is in principle an arbitrary vector. Quasinormal modes would come from imposing boundary 
conditions for the scalar field configuration at the horizon and at the boundary: We ask for incoming 
modes boundary conditions at r = r+ while we ask the field to vanish at r = oo. To implement these 
conditions it is useful to consider the following coordinate transformation 



y 



(21) 



for which the range r = [r+, +00) is mapped to y = [0, 1). This coordinate transformation induces a 
change of variable in the equation for the radial part of the scalar field R{r) ^ R (y), which now has 
to fulfill the following equation 

d^R (y) {2 + {D- 4) y) dR (y) P {ujH^ + {uH^ + k^rl) y^ - {2u;^l^ + k^l + rnVj) y) 



dy' 



■ + ^-r7^ .'"' r' +— ^ ':\. ' .0 ' '^^Riy) = 0, 

(22) 
where k^ = k ■ k. This equation can be transformed into a hyper geometric equation, whose solutions 
for iojl"^ /r'^ ^ Z read 



R (y) = y" (1 - yf [A^ 2F1 (a, 6, c, y) + A2 y"^" 2^1 (6 - c + 1, a - c + 1, 2 



c,y) 



(23) 



where A\ and A2 are integration constants, 2-^1 (y) stands for the hypergeometric function, and 



zo;. 



/3 



a 



/3 



L» + /D2T4^^^ 



2r\ ' '" 4 

where the arguments of the hypergeometric function are 
1 



4r^ 
1 



2 + \/d2T4^^) r\ - 2iiofi + ^J{D - 2f r| - 4/2 (/tV^. + ^2/4) 



4r2 



2 + VD^ + 4m2/2) rl - 2iLofi - ^ {D - 2f r\ - 4P {k'^rl + ^2/4) 



c= 1 



iujt 



(24) 
(25) 
(26) 



One can use properties of the hypergeometric functions to express the solution as a linear com- 
bination of functions that are regular in other points. We prefer to use the basis above in which the 
problem of recognizing the ingoing mode at the horizon is simpler. 

Imposing ingoing boundary condition at the horizon, which in these coordinates is located at y = 0, 
requires A2 to vanish, since around the horizon the scalar field behaves as 



J/— >0 y—^O 



(27) 



where g := /^/2r^, and where we have used the fact that 2F1 (^j b,c,0) = 1. Therefore, setting A2 = 
gives the desired behavior, namely ingoing modes at the horizon. One is then left with the following 



solution for the radial part 

R {y) = Aiy- (1 - y)f 2^1 (a, h, c, y) . (28) 

In order to impose a boundary condition at infinity, y — )• 1, it is convenient to use the Kummer's 
identities for the hypergeometric function and express 2-^1 (y) in terms of combinations of 2F1 (1 — y). 
After performing this transformation the function (j28p reads 

R {y) = A^y'^ (1 - yf [^i si^i (a, 6, a + 6 + 1 - c, 1 - y) + 6 (1 - 2/)'^"^' 2^1 (c - a, c - 6, 1 + c - a - 6, 1 - y) 

(29) 

with 

r(c)r(c-a-6) r(c)r(a + b-c) 

Ci = TT7 TTTT 77 and 6 = t- . .^,,. ■ (30) 

r (c — a) r (c — 0) r (a) r (o) 



Now, we impose Dirichlet boundary conditions at infinity: Expanding the expression (j29p when 
y — 7- 1, one has the following leading order contributions coming from each of the terms inside of the 
square bracket of (p9|) 

$ ~, ei (1 - y)^ +6(1- 2/)'^+^"'^"^ ~ 6r-^+ + e2r-^- , 

with 

A± := ^—^ . (31) 

Notice that this gives the solutions to equation ([3]). For w?' > the contribution given by r~ 

diverges at infinity, and then in order to impose Dirichlet boundary condition $ ~ 0, one is forced 

to impose 

_ T{c)T{a + b-c) _ 

^'- r(a)r(6) -°- (^^^ 

This is fulfilled provided a = —n or b = —n with n G Z>o with c — a — h ^ "L. These conditions give 
rise to the following spectrum 



^ Pk"^ + r\[D + m'^f + in {n + 1) + (2n + 1) VD'^ + Am'^P 

oj = --T !^ , ^ , (33) 

P VD^ + 4m2?2 + 2 (2n + 1) 

where we have reinserted r^ . In terms of the dimension A± this reads 

_ i Pk^ + rl (A_(A_ - 2(2n + 1) - D) + 2(2n + ■D)(n + 1)) 
'^~ Z3 2(2n + 1) + 1? - 2A_ ' ^ ' 



In the particular case D = A, the frequencies (j33p coincide with frequencies uji of reference [26]. 

At the special points c—a—b € Z, that is when ^/D"^ + AmPP /2 € Z, the solutions of hypergeometric 
equation develop logarithmic dependence. For m? > 0, the frequencies (j33p are purely imaginary and, 
besides, the imaginary part is strictly negative. This implies that there are no oscillatory modes and 
the field is purely damped. In this case the leading term in the asymptotic behavior of the scalar field 



is given by $ ~ r +, that we will call strong asymptotic behavior. 

J/-5.1 

Notice that for - {D/2f < m^f < both A± > 0, and so both branches r~^±, fulfill the Dirichlet 
boundary condition with a continuos spectrum. This is similar to the case of AdS spacetime [15] for 
fn\pl'^ < m^/^ < 0, where "m^Qp := —{D — 1)^/(2/)^ stands for the Breitenlohner-Preedman mass [25] . 
Notice also that the bound — {D/2) < m^t^ in the asymptotically Lifshitz black hole with z = 2, 
corresponds to the AdS-Breitenlohner-Preedman bound in dimension D + 1. 

3.2 Scalar field and Lifshitz soliton 

In three dimensions, the metric (|18p allows one to construct a soliton by performing the following 
Wick rotation both in the spatial and the time coordinates; namely 

il ~ ^ il^ , . 

X -^ — t and t -^ -^X . (35) 

The change in the radial coordinated r = r^ cosh p in the three dimensional version of (jlSp , supple- 
mented by the double Wick rotation (j35p . gives the soliton metric 



ds"^ = f [- cosh^ p dP + dp^ + cosh^ p sinh^ p dX"^] , (36) 

where p > 0, and the horizon of the black hole is mapped to the center of the everywhere regular 
soliton. Asymptotically, this metric has an anisotropic scaling symmetry with z = Zgouton = 1/2, and 
its Euclidean continuation is diffeomorphic to the Euclidean continuation of the black hole. Given the 
smoothness of the solution (j36p and its lack of integration constants, it would be natural to regard 
this spacetime as a groundstate. Some of the implications of this in the case z = 2> have been already 
explored in [28], where it was shown that considering as a background the soliton constructed out 
from the z = 3 black hole in three-dimensional massive gravity, it is possible to give account of the 
entropy of the black hole by means of a generalized version of Cardy's formula. 

It is therefore natural to raise the question about the stability of the solitons constructed in this 
manner. In particular, if one considers a scalar field perturbations on the soliton (j36p . it is possible to 
integrate the field configurations in terms of hyper geometric functions. Then, by imposing regularity 
at the origin [p = 0) and vanishing field at infinity {p — )• oo), one finds the normal modes of the field 
on the soliton. Nevertheless it is possible to avoid this computation since requiring ingoing boundary 
condition at the horizon of the z = 2 black hole (jlSp in three dimensions corresponds to requiring 
regularity of the scalar field at the origin of the soliton (j36p , and therefore the normal modes on the 
soliton can be obtained directly from the double Wick rotated quasinormal frequencies ()33p on the 
black hole in three dimensionqj. This is, performing the changes 

k -^ —iujsol and io -> -r^ikgoi , (37) 



The same occurs for BTZ where the sohton obtained after the Wick rotation is AdSa. 



for D = 3 in (j33p . and then solving for the frequencies on the soUton UgoU one gets the following 
spectrum 

Usoi = ± ((2n + |A;,o;| + 1) \/9 + AmH'^ + 3 + m'^f + 4n (n + \ksoi\ + 1) + 2\ksoi\) , (38) 



where n S Z>o, and fc^oj is the momentum along the X direction in (j36p . 

Note that the normal frequencies in (138p are real provided m^ > —9/(4/^), so the stability of the 
scalar field propagation on the soliton is guaranteed on this range. 

3.3 Scalar field stability for generic z 

We observe from (j33|) that the imaginary part of the frequency is actually negative, namely Im{uj) < 0. 
As we will show below, this result could have been shown by a more general argument, which comes 
from adapting an argument of [H], valid for asymptotically AdSu black holes, to generic values of 
z. More precisely, having Im{uj) < for Dirichlet boundary conditions is a generic feature of the 
asymptotically Lifshitz black holes of the form (J16p witlo F (r) = 1 — r'^/r'^, and for arbitrary values 
of the dynamical exponent z. To see this, let us consider the ansatz 

d^l = -^nr)dt' + ^a^^ + r'dJ:l_, , (39) 

with d'E'j^_2 being the metric on the {D — 2)-dimensional base manifold. By a simple change of 
variables, metric (j39p can be rewritten as follows 

dsl = -r^^F (r) dv^ + 2r'~^dvdr + r^dS^.g , (40) 

where we have fixed P = 1 for short. 

Now, consider again the scalar field equation (|19p and propose the following form for its solution 



$ = e-^'^''i? (r) y (S) , (41) 

where Y (S) is a harmonic function on S, i.e. V|y (S) = —k'^Y (S) where V^ is the Laplacian 
intrinsically defined on S/)_2j in our case T^d-2 '■^ M . Using this, one obtains 



y(s) 



dr dr \ dr ' 



0. 



2~ D 

If one further defines R (r) = r^^ip (r), this equation becomes 



d(f{r)di.\ 2^^^ _ ^(,) ^ ^ , (42) 



dr \r^ ^ dr I dr 



^The same argument works without assuming the special form F(r) = y = 1 — r\/r'^ but merely requiring F' [r) > 
outside the horizon. 



10 



where / (r) = r^^ ^ (r^ — l) , with r^ = 1, and where the effective potential is 

V{r) = -^ {{D -2){D-6 + 2z) r^^'^ (r^ _ l) + 4 ((D - 2) + m^) r^^ + Ak'^r'^'-'^) , (43) 

Notice that for z > 2 this potential is strictly positive outside the black hole. Notice also that in 
the case z = 1, this potential is positive for r > 1. 

By manipulating equation (I42p and integrating by parts one can show that the following equation 
holds 



where Dirichlet boundary conditions for the scalar field were imposed at infinity. Condition ([5]) 
guarantees the convergence of the left hand side of (|44p . Then, since the left hand side of (j44p is strictly 
positive, it demands Im (co) < 0, and then we conclude the stability of the scalar under perturbations 
respecting these boundary conditions, provided the generalized Breitenlohner-Freedman bound © 
is obeyed. This analysis includes, in particular, the z = 3 black hole solution of three-dimensional 
massive gravity, whose quasinormal modes were recently studied in |27j . 

4 Conclusions 

In this paper we have studied scalar field perturbations in asymptotically Lifshitz black holes with 
dynamical exponent z = 2 in arbitrary dimension D. These solutions appear, for instance, as exact 
solutions of higher-curvature theories of gravity, which consist of augmenting Einstein-Hilbert action 
with terms that are quadratic and cubic in the Riemann tensor. 

For these black holes, we have explicitly solved the quasinormal modes for scalar field perturbations 
that obey suitable boundary conditions. We have computed the frequencies of these perturbations 
explicitly, finding that they are purely imaginary, with negative imaginary part. This implies the 
stability of these geometries under scalar field perturbations that respect the mentioned boundary 
conditions. The stability of these Lifshitz black holes under scalar field perturbations, on the other 
hand, can be anticipated by a more general argument: By adapting a computation previously done 
for asymptotically AdS^j black holes, we gave a concise argument that proves the stability of a more 
general class of black holes with arbitrary value of z provided the generalized Breitenlohner-Preedman 
bound (BF) is satisfied. 

By performing double Wick rotations of the black hole (|18p . we have constructed an everywhere 
regular solitonic solution in three dimensions. This solution has an asymptotically anisotropic sym- 
metry. We have obtained the normal frequencies for the scalar field on this metric by performing a 
double Wick rotation of the quasinormal spectrum of the black hole. This is consistent since asking 
for an ingoing boundary condition for the scalar field on the horizon of the black hole corresponds to 
asking for regularity of the field at the center of the soliton. The spectrum obtained in this way shows 
that the scalar field oscillates without any damping neither exponential grow, provided the generalized 
BF bound is fulfilled. This analysis can be extended to higher dimensions, where in order to construct 



11 



a soliton one must single out one of the flat directions on the metric (jlSp . and then apply a double 
Wick rotation between the selected coordinate and time. Performing the corresponding Wick rota- 
tion on the quasinormal spectrum (j33p itself in arbitrary dimensions gives the normal modes for the 
scalar perturbation on the soliton with frequencies that are real provided the generalized BF bound 
is fulfilled. It would be interesting to see if the nonlinear instability of AdS under the collapse of a 
backreackting scalar field [5U], also appears in the case of the asymptotically Lifshitz backgrounds. 

Before concluding, we would like to emphasize that in this article we have only considered the 
quasinormal mode computation for scalar fields, and even when from the holographic point of view such 
computation is interesting by its own right, concluding stability of asymptotically Lifshitz spacetimes 
would also require to consider, in particular, spin-2 excitations. Being a higher-curvature theory, 
computing spin-2 perturbations of an action like (|17p about an asymptotically Lifshitz geometry is 
certainly cumbersome, and since the equations of motion are of fourth order, the results would in 
general exhibit ghosts. Still, particular higher-curvature models like the one introduced in [13] could 
likely result to be free of ghosts about Lifshitz black hole solutions as well. This is an interesting 
question that requires further study. 
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grants, from ANPCyT, CONICET and UBA, Argentina. The authors thank Pedro Alvarez, Eloy 
Ayon-Beato, Francisco Correa, and Bertha Cuadros-Melgar for comments. G.G. and M.L. thank the 
hospitality of Universidad Austral de Chile during their stay. After our paper appeared in arXives, 
we were informed by Oscar I. Fuentealba Murillo and his collaborators that, in an unpublished paper, 
he and collaborators performed computations similar to those of our section 3.2. 

5 Appendix 

The coupling constants of action (|17p for the theory to admit z = 2 Lifshitz black holes (J16p in 
arbitrary dimension D are given by 

A = -'^—-D{^D^ - 26D^ + W7D^ - 20D^ - SOD"^ + 368D + 208), (45) 

3Cd 

bi = ^{3D^-15D^ + 6D^ + 10iD + 32), (46) 

^ {D-1){D^ -0^^-200^ + 108D'^ + 208D + 6A), (47) 



Cd {D + 2) 
"^^ = "^^^^ + ^^^^ - 8)(D - l)\ (48) 

C2 = -:t- {D^ - ISL'^ + 8D + 4) , (49) 

8cr/^ , , 

C3 = ^^i4D + l), (50) 
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where Cd ■= D^ - 8D^ + 531?'^ - 130D^ + 124L>2 + 296D + 64, which has no integer roots. 
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